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TOWARDS PRECISION DETERMINATION OF UNINTEGRATED 
PARTON DENSITY FUNCTIONS 
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First attempts are described to determine the unintcgratcd Parton Density Function of the gluon 
from a fit to measurements of the structure function F2(x,Q 2 ) and also ^(ajjQ 2 ) measured at 
HERA. Reasonable descriptions of both structure functions are obtained, but the gluon densities 
determined are different. 
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1. Introduction 

Unintcgratcd parton density 

functions (uPDFs) are best suited to study 
details of the hadronic final state in high en- 
ergy ep and also in pp collisions (for a review 
see [1~7]). In general, the production cross 
section for jets, heavy quarks or gauge bosons 
can be written as a convolution of the uPDF 
A{x,k\,q) with the partonic off-shell cross 
section a(xi, k±), with Xi,k± being the lon- 
gitudinal momentum fraction and the trans- 
verse momentum of the interacting parton i 
and q being the factorization scale. For ex- 
ample the cross section for ep — > jets + A" can 
be written as: 



dErdr] 



E 



dxi dQ d . 



[dk\x t A{xi, kj_,q)] a(xi,kj_) 



At high energies, the gluon density is dom- 
inating for many processes, therefore here 
only the gluon uPDF is considered. It has al- 
ready been shown in [8] , that the predictions 



of the total cross section as well as differen- 
tial distributions for heavy quark production 
at HERA and the LHC agree well in general 
with those coming from fixed NLO calcula- 
tions. However, the details depend crucially 
on a precise knowledge of the uPDF. There- 
fore precision fits to inclusive and exclusive 
measurements have to be performed to de- 
termine precisely the free parameters of the 
uPDF: the starting distribution function at 
a low scale qo ~ 1 GeV as well as parameters 
connected with a s and details of the splitting 
functions for the perturbative evolution. 

The previously available uPDFs (a 
overview is given in [3-5]) were only adjusted 
to describe F2, but no real fit of the pa- 
rameters of the starting distribution was per- 
formed, neither the experimental uncertain- 
ties were treated as is done in global QCD 
analyses in the collinear approach (for exam- 
ple from CTEQor MRST). 

Here attempts to determine the uPDF of 
the gluon from pQCD fits using the CCFM 
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evolution equation [9] . to the inclusive struc- 
ture function F2(x, Q 2 ) as well as F£(x, Q 2 ), 
as measured by the HERA experiments, are 
described. 

2. The method 

The unintegrated gluon density is de- 
termined by a convolution of the non- 
perturbative starting distribution Ao (%) an d 
the CCFM evolution denoted by A (x, k±, qj: 

xA(x,k±,q) = / dx'Aoix') ■ — 
J x' 

In the perturbative evolution the gluon 
splitting function P gg including non-singular 
terms (as described in detail in [10]) is ap- 
plied. 

The distribution Ao is parameterized at 
the starting scale q~o by: 

xAo(x) = Nx~ B ' ■ (1 - xfo (1 - D g x) 

The parameters N g , B g ,C g , D g of Ao are de- 
termined in the fit, which is based on the MI- 
NUIT package [11] and the extension by [12], 
such to minimize the \ 2 defined by: 

^ / t I g2 stat _j_ g2 uncor J 

with T being the theory value and D the 
measurement with the corresponding sta- 
tistical and uncorrelated systematic uncer- 
tainty. 

3. The structure function 
F 2 (x,Q 2 ) 

The measurement of F2(x, Q 2 ) [13] is used in 
the range x < 0.005 and Q 2 > 5 GeV 2 to 
determine the uPDF. The parameters of Ao 
were investigated separately, and it was con- 
firmed, that the data do not constrain the 
parameters C' g and D g , which were there- 
fore fixed to C g = 4 and D g = 0. The 
starting distribution Ao was parameterized 



at go — 1-2 GeV. The running coupling 
Q! s (/i) was used in the 1-loop approximation 
in the region (j, > qo, and was kept fixed 
at a s (qo) for fi < q~Q. A acceptable fit to 
the measured Fi values was obtained with 
x 2 /ndf = 111.8/61 = 1.83 using only sta- 
tistical and uncorrelated systematical uncer- 
tainties (compare to x 2 / n( V ~ 1-5 in the 
collinear approach at NLO). In Fig. 1 the 
measurement is compared to the prediction 
of the structure function F2(x,Q 2 ) as ob- 
tained from the fit. Since in the cross section 
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Fig. 1. The structure function i*2(^, Q 2 ) as a func- 
tion of Q 2 for different values of x. The data points 
are from [13]. The shaded area shows the region 
which is not used in the fit. 

a product of the gluon density and a s enters 
(as well as a s enters in the pQCD evolution) 
also the dependence on the choice of A qc d was 
investigated (Fig. 2). A clear preference for 
^qcd ~ 0.13 is observed, which corresponds 
to a s (mz ) = 0.118 in the 1-loop approxima- 
tion. 
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Fig. 2. x 2 / n df dependence of the F2 fit on the 
choice of A gcl j 

4. The structure function 
F£(x,Q 2 ) and the uPDF 

The structure function F^x^Q 2 ) is directly 
sensitive to the gluon density. The measure- 
ments of [14] were used to determine the 
uPDF in the range Q 2 > 1.5 GeV 2 . A ac- 
ceptable fit to the measured F-f values was 
obtained with x 2 /ndf = 18.8/20 = 0.94 
using statistical and systematic uncertain- 
ties. In Fig. 3 the measurement is compared 
to the prediction of the structure function 
F^x^Q 2 ) as obtained from the fit. 

Since the fits to F 2 (x, Q 2 ) and F 2 c (ie, Q 2 ) 
were performed independently, it is interest- 
ing to compare the resulting unintegrated 
gluon distributions (Fig. 4). The small x be- 
havior of the gluon distribution is very differ- 
ent. It has been explictly checked, that the 
measured F2 (x, Q 2 ) cannot be described with 
the uPDF obtained from the fit to F£(x, Q 2 ) 
and vice versa. 

5. The structure function 

Using the different uPDFs described in the 
previous sections we calculate the longitu- 



dinal structure function Fl(x,Q 2 ) at fixed 
W in the framework of [16]. The predic- 
tions are compared with the measurements 
at fixed W of HI [15] in Fig. 5. The general 
trend is nicely described, however the pre- 
dicted shape of Fl at medium Q 2 is different 
and precise measurements would be very wel- 
come. 
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Fig. 4. The unintcgratcd gluon density as obtained 
from the fits to F2(x, Q 2 ) and F£{x, Q 2 ) as a function 
of x for different values of kj_. The gluon density is 
evolved to a scale q = 10 GeV. 
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Fig. 5. The longitudinal structure function 
Fl(x, Q 2 ) at fixed W as measured by [15]. The solid 
curve is solid curve is the prediction using the fit to 
i*2 whereas the dashed curve uses the uPDF obtained 
from a fit to F§ ■ 
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